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We discuss the current status of the resonant spin-flavor precession (RSFP) solution to the solar 
neutrino problem. We perform a fit to all the latest solar neutrino data for various assumed magnetic 
field profiles in the Sun. We show that the RSFP can account for all the solar neutrino experiments, 
giving as good fit as other alternative solutions such as MSW or vacuum oscillation (Just So), and 
therefore can be a viable solution to the solar neutrino problem 
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I. INTRODUCTION 
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The solar neutrino anomaly has been more and more strongly established both in its experimental ^-[^] as well as 
in its theoretical aspect. In fact, both have presented a very dynamical evolution. From one side, theoretical 

predictions, i.e., the standard solar models (SSM), have been refined by including several mechanisms such as helium 
qq ' diffusion and an updating analysis of the Sn astrophysical factor |^,[l(|. One can see in ref. jllj that theoretical 
predictions obtained by different SSMs, which are developed independently, are in good agreement with each other. 
t^J" ' Moreover, it has been shown that the SSM is in excellent agreement with the helioseismology [||. 

On the other hand, experimental data have become more accurate due to the calibration of the experiments, 
more statistics and the existence of the new generation SupcrKamiokande experiment and its solar neutrino spectral 
observations [El . Consequently, final numbers related to the solar neutrino deficit have also evolved significantly. The 
most updated experimental data as well as theoretical predictions are shown in Table I. One can show that these 
observed solar neutrino data are in strong disagreement with the ones predicted by the SSM |0,Q. Moreover, this 
Qh, conclusion does not depend on any details of the SSM. 

Solutions to the solar neutrino problem rely on different phenomena jl3| which deplete the number of observable 
neutrinos at the Earth: neutrino oscillations in vacuum JlJ], resonantly enhanced matter oscillations [flEf , resonant 
' spin-flavor precession phenomenon (RSFP) |l(],[l7j and flavor mixing by new interactions Jl^Jl9|| . The capability of 
each one of these processes to make compatible solar neutrino predictions and observations and therefore to find a 
solution to the solar neutrino anomaly have been updated from time to time pp[ . In this paper we investigate the 
current status of the RSFP scenario [£lj . We believe that it is worthwhile to reanalyse this mechanism in the light of 
new solar neutrino data as well as the new SSM. We also discuss briefly how the solar neutrino spectral observations 
in SupcrKamiokande are affected by the RSFP mechanism. 



TABLE I. Observed solar neutrino event rates used in this analysis and corresponding predictions from the reference standard 



solar model M. The quoted errors are at la. 



Experiment Data ±(stat.) ±(syst.) Ref. Theory 




Units 


Homestake 2.56 ± 0.16 ± 0.15 

SAGE 69.9± s 7 ° T ±li 
GALLEX 76.4 ± 6.3+* j 
SuperKamiokande 2.44 ± +0.05±^6 


1 
' 

1 
5 


7.7±i-% SNU 
129±| SNU 
129+^ SNU 
5.15+°' 9 72 10 6 cm-V 1 
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RSFP mechanism is very sensitive to the magnetic profile in the Sun and, in fact, several possible scenarios for the 
magnetic strength have been proposed by different authors J2^-|2(|. We consider several possibilities which include 
the main qualitative aspects of the magnetic profile in the Sun previously invoked as a solution to the solar neutrino 
anomaly. Using the minimum \ 2 method to compare theoretical predictions of the RSFP phenomenon and the 
experimental observations we conclude that very good fits can be obtained for the average solar neutrino suppression, 
if intense magnetic fields in the solar convective zone are considered. 



II. RSFP MECHANISM 



Assuming a nonvanishing transition magnetic moment of neutrinos, active solar neutrinos interacting with the 
magnetic field in the Sun can be spin-flavor converted into sterile nonelectron neutrinos p7| , p8[ (if we are dealing with 
Dirac particles) or into active nonelectron antineutrinos |2!|] (if the involved particles are Majorana). In both cases the 
resulting particles interact with solar neutrino detectors significantly less than the original active electron neutrinos 
in such a way that this phenomenon can induce a depletion in the detectable solar neutrino flux. 

Spin- flavor precession of neutrinos can be resonantly enhanced in matter (T^jlTj] , in close analogy with the MSW 
effect [ jl5| . In this case the precession strongly depends on the neutrino energy and provokes different suppressions 
for each portion of the solar neutrino energy spectrum. Therefore RSFP provides a satisfactory description [^l] ^Oj 
of the actual experimental panorama : all experiments detect less than the theoretically predicted solar neutrino 
fluxes Q and different suppressions are observed in each experiment, suggesting that the mechanism to conciliate 
theoretical predictions and observations has to differentiate the different parts of the solar neutrino spectrum. 

For simplicity, we consider two generation of neutrinos, electron neutrino and, for e.g., muon neutrino (which could 
be replaced by tau neutrino in our discussion). Furthermore, we assume that the vacuum mixing angle is zero or small 
enough to be neglected. (See ref. ^] for the case where RSFP and flavor mixing simultaneously exists.) The time 
evolution of neutrinos interacting with a magnetic field B through a nonvanishing neutrino magnetic moment \i v in 
matter is governed by a Schrodinger-like equation 



dt 



ii R Am 2 



(1) 



where v e and Vu R are active electron neutrinos and muon antineutrinos, respectively, Am 2 



is their 



squared mass difference and E is the neutrino energy, a Ve = Gf(2A^ 6 — N n )/^/2 and a„ — Gf JV„/\/2, with N e and 
N n being electron and neutron number densities, respectively. In eq. (0) we are assuming that neutrinos are Majorana 
particles. For the Dirac case, the spin-flavor precession involves v e <-> v s , where v s is a sterile neutrino and a Vs = 0. 
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FIG. 1. Matter potentials as a function of radial distance from the solar center are plotted. The solid and dotted curves 
correspond to the Majorana and Dirac case, respectively. The dashed line correspond to fi„B = 10 _11 /is x 100 kG. 



III. ANALYSIS 

In order to obtain the survival probability we first numerically integrate the evolution equations (0) varying matter 
density in the Sun Q for some assumed profiles of the magnetic field which will be described below. Next, using 
the solar neutrino flux in ref. ||, we compute the expected solar neutrino event rate in each experiment, taking into 
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account the relevant absorption cross sections |J for 71 Ga and 37 C1 experiments as well as the scattering cross sections 
for v e -e~ and i^-e~ reactions including also the efficiency function for the SupcrKamiokande experiment in the same 
way as in ref. M]. We note that in this analysis we always adopt the solar model in ref. § as a reference SSM. 
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FIG. 2. Various magnetic field profiles used in this work. For each field (B) is defined as the average of the field over the 
region where B(r) is not zero. 

A. Assumptions 



We do not take into account the production point distribution of neutrinos. This can be justified by two main 
reasons. For the relevant values of Am 2 our analysis shows that (see below) the resonance position lyes always outside 
the neutrino production region (r/i?suN < 0.3) because the solutions we find implies Am 2 < 10 -7 eV 2 . In order to 
see this we plot in Fig. 1 the matter potential as a function of the radial distance from the solar center. 

If neutrinos are considered Majorana particles, V ma tter = cii/ e — o-v^ , and if they are Dirac, Vmatter = o-u e ■ It is shown, 
in the right ordinate in the same plot, also the value of the quantity Am 2 to which a resonance is found. For example, 
if Am 2 = 10~ 7 eV 2 (for E= 1 MeV), a resonance is localized in r/i?suN w 0.52. Therefore, it can be seen that the 
resonance position is also larger than r/i?suN = 0.3, outside the production region. Also, the range of [i v B we consider 
here is much smaller than the matter potential, a Vc and a„ M , in the production region. Again this fact is shown in 
Fig. 1, where the quantity fi„B = 10~ n /iB x 100 kG is presented. Therefore, from eq. ([!]) we observe that neutrino 
spin-flavor precession is prevented before neutrino gets the resonance region and the final survival probability is not 
affected by the production position. 

It is obvious from the evolution equations (|l|) that the RSFP mechanism crucially depends on the solar magnetic 
field profile along the neutrino trajectory. In the present paper we choose several different profiles which, we believe, 
cover in general all the previously |22|-[2q| analysed magnetic profiles which led to a solution to the solar neutrino 
anomaly. In Fig. 2, these magnetic fields are presented in their general aspects. The constant magnetic profile B\{r) 



jted in references [E4|, while the general aspects of the profiles B^(r) and B±(r) have already appeared in 
25|,|2^], and |2(|, respectively. We note that close to the solar surface (r > 0.95-R©) we have switched off the 
magnetic field for all the profiles we considered in this work. 
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B. Definition of \ 2 

The relevant free parameters in the RSFP mechanism are Am 2 and \i v B. Using the solar neutrino data given in 
Table I, we look for the region in the Am 2 — fJ, v (B) parameter space, where (B) denotes the average field strength 
defined as in Fig. 2, which leads to a solution to the solar neutrino anomaly by means of the minimum x 2 method. 
In this analysis we will use only the SuperKamiokande data Q without including the Kamiokande data [0 because of 
the larger statistics and the smaller systematic error in the SuperKamiokande experiment. We also note that we will 
use the combined value of the two 71 Ga experiments, 72.3 ± 5.6 SNU. 

Our x 2 is defined as follows, 

x 2 = - Rf'M^T'i*? - R T S )' ( 2 ) 

hj 

where run through three experiments, i.e., 71 Ga , 37 C1 and SuperKamiokande, and the total error matrix af^tot) 
and the expected event rates Ri are computed as follows. We essentially follow ref. |33] for the derivation of the error 
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matrix and to describe the correlations of errors we used in this work. The expected experimental rates in the absence 
of neutrino conversion is given by, 



Rt = c vh (* = Ga ' C1 > SK )' ( 3 ) 

3 

where Cy is the cross section coefficients and 4>j is the solar neutrino flux. In Sec. IV where we consider the case with 
neutrino conversion we use the coefficients CV,- determined by properly convoluting the conversion probability in the 
integration over the neutrino energy spectrum for each detector and neutrino flux. In this work we consider neutrino 
from pp, pep, 7 Be, 8 B, 13 N and 15 sources and neglect other minor flux such as 17 F and hep neutrinos. 
The total error matrix 4 is the sum of the theoretical afj(th) and experimental one cr^-(exp), 

4 (tot) =4(th) + 4(exp). (4) 

The theoretical error matrix can be further divided into the one coming from the uncertainties in the cross sections, 
afj (cross) and the one coming from uncertainties in the solar neutrino flux, 4 (flux) , 

4 (th) = 4 (cross) + 4 (flux) . (5) 
The cross section error matrix 4 (cross) can be calculated by, 

4(cross)^, ± ^.^. Al^AlnC, 



k,i=i 



= %^(i? lfc AlnC7 jfc ) 2 , (6) 



fe=i 

where R ik = C ik 4>k- 

1 



On the other hand, the flux error matrix 4 (flux) can be calculated by, 

k,l=l m=l 
6 11 

= RikRji ^(Aln<fe) m (Aln^)m, (7) 

fc,i=l n»=l 

where (Aln^) m is the fractional uncertainty of the k-th neutrino flux coming from the uncertainty in the m-th input 
parameter (S11, 533, S34, Sn, Z/X, opacities, 151,14, luminosity, age, diffusion or 7 Bc + e~ capture rate) which are 



obtained by the computer code exportrates.f, available at URL http://www.sns.ias.edu/~jnb/. We note that in eq. 
(fjj) the product (Aln^/ c ) m (Aln0;) m takes positive (negative) value when fc-th and Z-th fluxes are positively (negatively) 
correlated to each other with respect to the variation of the m-th input parameter. 
The experimental error matrix is given by, 

4( ex P) = hWji ( 8 ) 

where crjj (i,j — Ga, CI, SK) stands for the combined error in each experiment. 
In Table II we show the correlation matrix pij defined as, 

4( tot ) , Q , 

Pij = , (9) 
^4 (tot) a*, (tot) 

TABLE II. The correlation matrix pij obtained from eqs. (M-(fe>|). 



Experiment 


Correlation matrix 






Ga 


1.00 






CI 


0.497 


1.00 




Super-Kam 


0.486 


0.952 


1.00 
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IV. RESULTS 



Now we compute the spin-flavor conversion probability, by numerically integrating the evolution eq. (1), assuming 
fi v = 10 _11 /is as a reference value, which is slightly below the present experimental upper bound |35| ]. Hereafter we 
always assume this value of magnetic moment |m| in this work but as it is clear from eq. (1) if /i„ is assumed to 
be smaller by certain amount, the same effect can still be obtained by simply increasing the value of magnetic field 
strength properly so that the product [i v B would not be changed. 

In Fig. 3 (a), (b), (c) and (d) we show the contour plots of the survival probability P(v e — > v e ) = \{v e (t)\v e (0))\ 
in the Am 2 — (B) parameter space, for the magnetic profiles we sketched in Fig. 2, B\(r), B 2 {r), B 3 (r) and B±{r), 
respectively. We note that the field B% gives very different probability contours from the other profiles, which will be 
also reflected in the final allowed region (see below). 

Including now the experimental observations on the solar neutrino signal above shown in Table I, we can determine 
the region in the Am 2 — (B) parameter space which leads to a RSFP solution to the solar neutrino problem for a 
specified confidence level. We present the Am 2 — (B) parameter region which can account for all the solar neutrino 
data, at 90, 95 and 99 % C. L. in Figs. 4 (a), (b), (c) and (d), for the magnetic profiles B\{r), B 2 (r), B 3 (r) and B±{r), 
respectively. 

We observe from Figs. 4 (a) to (d) that a solution to the solar neutrino problem can be found when (B) > few times 
10 kG and Am 2 is the order of 10~ 8 to 10~ 7 eV 2 for any of the magnetic profiles used in this work. Nevertheless, the 
quality of the fit, measured by the minimum x 2 criterion, varies a lot. The poorest fit is obtained when the continuously 
decaying magnetic field profile B^ is used, with Xmin = 6-1 for one (three data points - two free parameters) degrees 
of freedom. Better fits are obtained when the B\ (uniform) and B 2 (large triangle) fields are employed showing 
Xmin — 2-0 and 1.8, respectively. And the best fit appears when the triangular field in the solar convective zone, B$, 
is employed, with a rather small value x m in = 0.13. For this profile, we note that, as expected from Fig. 3 (c) we 
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have several local best fit points also indicated in Fig. 4 (c) by the open circles whose corresponding Xmin are ' f rorn 
left to right, 2.3, 0.29 and 0.19. 
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FIG. 4. The Allowed RSFP solution to the solar neutrino problem. The parameter region allowed at 90, 95 and 99 % C. 
L. are shown in (a), (b), (c) and (d) for the magnetic field profiles, B\, B2, -B3 and B4, respectively, sketched in Fig. 1. We 
indicate best fit points by filled circles. In (c) we also indicate, by the open circles, the local best fit points inside each island 
delimited by 90 % C.L. curves. 



The reason why we are getting very good fit for B3 is that this profile can provide the required suppression patterns 
of various neutrino flux implied by latest the data Jl2| as discussed in ref. |^5fl . First we note that low energy pp 
neutrinos are not so suppressed because the resonance positions are located in the inner region where the magnetic 
field is zero or small (see Figs. 1 and 2). However, intermediate energy 7 Be neutrinos can be strongly suppressed due 
to the rapid increase of the field at the bottom of the convective zone since their resonance position is in a slightly 
outer region than the pp one. On the other hand, high energy 8 B neutrinos are moderately suppressed because their 
resonance positions are closer to the solar surface than the 7 Be ones, where the field is decreasing. The best fitted 
values of (B) and Am 2 as well as Xmin obtained from these different profiles are summarized in Table III. 

TABLE III. The best fitted parameters and Xmin f° r the Majorana case. Dirac case is presented in the parentheses. 



Profile 


(B) (kG) 


Am 2 (10~ 8 eV 2 ) 




1 


50.6 (40.9) 


3.5 (2.3) 


2.0 (6.2) 


2 


47.1 (40.8) 


6.1 (4.6) 


1.8 (5.7) 


3 


118 (69.4) 


1.5 (1.2) 


0.13 (1.3) 


4 


82.9 (81.6) 


8.1 (6.6) 


6.1 (11.4) 
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We have repeated the same analysis also for the Dirac neutrino case. We, however, do not show the plots for the 
allowed region since they are rather similar to what have been presented above, if the same magnetic field profile 
is assumed. Instead, for the case of Dirac neutrinos, we only present the best fitted parameters and Xmin m ^ ne 
parentheses in Table III. 

We see from this table that, the Dirac case always leads to a worse fit if the same magnetic field profile is assumed. 
To understand this we should note that for the Dirac case, iVs are converted into the right handed muon (or tau) 
neutrino, which do not contribute to any of the solar experiments including the water Chrenkov experiment p7| . 
In contrast in the Majorana case, converted right handed neutrino D^s do contribute to the signal observed in the 
SupcrKamiokande detector. This makes it difficult to conciliate the difference between the SuperKamiokande and 
37 C1 data. 

Let us now comment about the possibility of having such strong magnetic field in the Sun. While there is no generally 
accepted theory of solar magnetic field, it is possible to bound the field strength from very general arguments. It can 
be shown that the magnetic field less than 10 6 kG in the solar core or less than 10 4 kG in the solar convective zone, 
will hardly affect the thermal structure and nuclear reaction processes well described by the standard solar model. 
These values come from the requirement that the magnetic pressure should be much smaller than the gas pressure, 
and can be regarded as the most generous upper limits of the magnetic field inside the Sun. More stringent bounds 
on the magnetic field in the convective zone are found in refs. []35| , j36| where the discussion is based on the non- linear 
effects which eventually prevents the growth of magnetic fields created by the dynamo process. Naive limit can be 
obtained by estimating the required field tension necessary to prevent a fluid element from sinking into a magnetically 
stratified region, so that the magnetic flux would not be further amplified. By equating the magnetic tension to the 
energy excess of a sinking element at the bottom of the convective zone, Schmitt and Rosner |35j] obtained ~ 10 kG 
as an upper bound for the magnetic field, which is of the order of the magnitude we need to have a good fit to the 
solar data by RSFP mechanism for the reference value of magnetic moment, /i„ = 10~ n fiB- 

Finally, we briefly discuss how the recoil electron energy spectra in the SuperKamiokande detector will be affected 
by the RSFP mechanism (38). In Fig. 5 (a) we plot the electron neutrino survival probabilities as a function of neutrino 
energy using the best fit parameters. In Fig. 5 (b) we plot the recoil electron energy spectra divided by the standard 
prediction expected to be observed in the SuperKamiokande detector, using also the best fit parameters as in Fig. 5 
(a). In Fig. 5 (b) we also plot the latest data from SuperKamiokande ||. As we can see from the plot the observed 
data indicate some distortion mainly due to the last three data points in the higher energy bins. We, however, note 
from this plot that it seems difficult to exclude, at this moment, any of our predicted spectra expected from different 
field profiles, because of the experimental errors. We have to wait for more statistics and more careful analysis from 
the experimental group before drawing any definite conclusion. 
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FIG. 5. We plot in (a) electron neutrino survival probability as a function of energy with the best fitted parameters for 
various field profiles. In (b) we plot recoil electron energy spectra expected from RSFP scenario using our best fit parameters, 
divided by the SSM prediction. The SuperKamiokande data are also shown by the filled circles with error bars. The last data 
point includes the contribution from the electrons with energy larger than 14 MeV. 



V. CONCLUSIONS 



We have reanalysed the RSFP mechanism as a solution to the solar neutrino problem in the light of the latest 
experimental data as well as the theoretical predictions. 
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We found that the quality of the RSFP solution to the solar neutrino anomaly crucially depends on the solar 
magnetic field configuration along the neutrino trajectory inside the Sun. We found that the best fit to the observed 
solar neutrino data, which seems to be even better than the usual MSW solution as far as the total rates are concerned, 
is obtained if intensive magnetic fiels in the convective zone is assumed, in agreement with the conclusion found in ref. 
p5[ , whereas the linearly decaying magnetic field gives the worst fit. We, however, note that the required magnitude 
of the free parameters involved in the process, i.e., the magnetic field strength multiplied by the neutrino magnetic 
moment n v (B) and the squared mass difference Am 2 , points to the same order, Hv(B) « few times 10 -11 /is • 10 kG 
and Am 2 ps few times 10 -8 eV 2 , for any of the field profiles assumed in this work. 

Our ignorance about the profile as well as the magnitude of the solar magnetic field makes this approach to the 
solar neutrino observation less predictive than its alternative approaches [ p^|Jl^Jl^| . Nevertheless the presence of this 
mechanism opens some interesting possibilities. 

One possibility is to look for any time variation of the solar neutrino signal [ |39| which can not be expected in 
other alternative solutions found in refs. [p^| , p^|Jl^ 1 . Any time variation of the solar neutrino signal which can be 
attributed to some time variation of the solar magnetic field can be a good signature of this mechanism. Although 
SuperKamiokande has not yet confirmed any significant time variation up to experimental uncertainty this possibility 
remains. 

Another possibility is to look for the solar T> e flux, which can not be produced in the usual MSW or vacuum oscillation 
case but can be produced in RSFP mechanism if the flavor mixing is included, produced by RSFP mechanism can 
be converted into v e by the usual vacuum oscillation. Ref. jh)| suggests to observe (or to put upper bound of) v e flux 
in the SuperKamiokande whereas ref. pH suggests to use low energy solar neutrino experiment such as Borexino or 
Hellaz. 

We finally stress that RSFP mechanism can still provide a good solution to the solar neutrino problem, comparable 
in quality to MSW or Just So solution, and is not excluded by the present solar neutrino data. 
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